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study the case when theses functions are polynomials. We construct 
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1 Introduction. 



Nonlinear ordinary differential equations appear in many branches of applied 
mathematics, physics and, in general, in applied sciences. 

By definition a real autonomous differential system is a differential system 
of the form 

x = v(s), xGM^ 

where the dependent variables real , the independent 

variable (time t) is real and functions v(a;) = (f ^(x), f ^(x)) are continuous 
functions in D C M^. 

Definition 1.1 The smooth function g and the relation g{x) = are said 
partial integral and invariant relation of the vector field v(x) respectively if 

dg{v)\g(^^)=o = 0. 



In this paper we are mainly interested in to study the differential system 
which possess a given set of invariant relations. 

It is always helpful to look at this problem from another point of view. 
In this paper, we take an alternative viewpoint of starting with a given set 
of invariant relations and determining the form of the system which has such 
a set as invariant set. [Hlj, ll2j,^, ilSj , [2] , [9] , [S] , [H] , [25] , ^ , |2T], 

m] 

This approach was first developed by Eruguin in the paper " Construction 
of the whole set of ordinary differential equations with a given integral curve" 
published in 1952 p!0]. In that article the author stated and solved the 
problem of constructing a planar vector field for which the given curve is its 
invariant. It is important to observe that Eruguin considered only one curve, 
moreover he didn't require that this curve was necessarily algebraic. 

Eruguin proved that the most general planar vector field v for which the 
given curve 

g{x,y) = (1.1) 
is its invariant curve generates the following differential equations 

( x = iy{x, y){g, x} + a{x, y) ,^ 
\ y = J^{x,y){g, y} + b{x,y) 

where and p, a, b are functions which we determine from the condition: 
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dg{w) = ^{x,y), $|g=o = (1.3) 

where 

/} = dxgdyf - dygd^f 

These Eruguin ideas were apphed in different areas. In particular Zubov 
in [27j| constructed the planar system with a given region of stability. 
Zubov constructed the following vector field 

X = fl{gu x} + gi{i{f, x} + giipdi) 

y = fiidi, y] + 9i{i{f, y] + 9i^d2) 

where /, 7, di, d2, ^ are arbitrary functions which he choose in such a away 
that 

di{gi{f, y] - f{gi, y}) + d2{gi{f, x} - f{gi, x}) = 1 

Under this condition Zubov proved that the following relations holds 

dgi{v) = giilif, gi}) + gj(p{di{gi, y} - d2{gu x}) 

rfG'i(v) = (^Gi, InGi^ — 

91 

Galliulin in |^12j determines the most general vector field in for which 
the given relations 

^7,(^,x^x^...,x^) = o, j = i,2,..,s< N 

are the invariant relations, where gi, g2, .., gs are smooth independent func- 
tions. The constructed system is the following 

1 ^ 

i = P XI ^^^^'^^ " grad^i + Y. 

where x = co/(x^, x^, x^), Y is an arbitrary vector orthogonal to the 
vectors 

d 

grad^ij- = col{digj, d2gj, dNgj), j = 1, S, dk = 

r is the Grama determinant, Tij are the minors of F and $1, are arbi- 
trary functions: 

$,l5,=o = 0, J = 1,2,..,^. 

The aim of this paper is to extend the Eruguin- Galliulin ideas to the case 
when the number of the given invariant relations is bigger than — 1. The 
results which we expose have been systematically developed in 
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2 Definitions and statement of the main re- 
sults 

In this section we constructed the most general stationary differential system 
from the given set of partial integrals. 

First of all we introduce the following concept and notations which we 
shall use below. 

Definition 2.1 

We call the vector field: 



dgi{di) 
dg2{di) 



dgi{d2] 
dg2{d2] 



dgnidi) dgM{d2) 

dgM +i{di) dgM +i{d2) 

dg^idi) dgN{d2) 

di 82 



dgi{dN) $1 
dg2{dN) $2 

dgM {On) $m 
• • dgM+i{dN) Am+1 

dgN{dN) Atv 





(2.1) 



the Eruguin-Galliulin Vector Fields, where gi, g2, ■■■,gN are smooth func- 
tions, $1, $2) '^'Af are the Eruguin functions and Xm+i, ^m+2, are 
arbitrary functions: 



k\gk=0 



0, k = l,..,M, 



(2.2) 

dgj{Y) = Xj, j = M + l,..,N, 

dgi, dg2, dgM are given independents 1-forms and dgM+i, dgM+2, dg^ 
are arbitrary 1-forms which we choose in such a way that 



T 



dgi{di) dgi{d2) ... dgi^d^) 
dg2{di) dg2{d2) ... dg2{dN) 



{^1, ^2, •••,^iv} 7^ (2.3) 



dgN^di) dgN{d2) ■ ■ ■ dgN{dN) 
The functions $1, $2, ^m we call the Eruguin functions [T2] . 
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We can identify the vector field fl2.ip with the first order differential sys- 
tem 

x = TM-^w, (2.4) 
where M and w are the matrices: 

f ({dgj{dk))j,k=i,2,..,N), 

y w = co/($i, Aaz+i, Aat). 

It is easy to show that the system (12. 4p admits the equivalent represen- 
tation 

= (^i{x\g2, ..gM+i--,gN} + ■■■ + ^A/lfl-i, ••, gu+u ••, S'a^} + ^^ (2.5) 
where 

= \M+i{gi--,gM, x\ gM+2--,gN}--- + XN{gi--,gM,gM+i--,gN-i, x^}, 

j = l,2,..,N. 

Clearly, the vector Y = col{Y^, Y^, ...,Y^) is orthogonal to the vectors 
gradgj, j = 1, 2, .., M, hence we obtain the Galliulin result |12] . 
Example 2.1 

We shall construct the Eruguin-Galliulin vector field for the case when 
the arbitrary functions 

gM+i, gM+2, ■■■,gM+K, n = m + k 

are such that 

dg{v)=Lg, 
dgu+iiy) = LgM+i + Lig 
i dgM+2{'v) = LgM+2 + ^ifl-M+i + ^25^ (2.6) 

dgM+Kiy) = LgM+K + LigM+K-i + ••• + Lxg 
where Li, L2, Lk, L are arbitrary functions and 

M 
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By introducing the functions Gi, G2, Gk '■ 

QM+j^GjQ, j^l,2,..,K, 

we obtain 

dGi{v) = Li 
dG2{v) = LiGi + L2 

_ dGK{^) = L^Gk-i + ... + Lk 

Clearly, the arbitrary functions Xm+i, Xm+2, Ajv = Am+k in this case we 
determine as follow 

Am+1 = g{LGi + Li) 
Am+2 = g{LG2 + LiGi + L2) 

^M+K — g{LGK + LiGk-1 + ■■■ + Lk) 

Let us introduce the 1-forms uji, a;2, ■■■iOJk '■ 

dGi — uji 
dG2 = GiUJi + UJ2 



I dGx — Gk-iOOi + ■■■ + ouk- 
After some straightforward calculations we prove that 

ujj = dTj, 



Hence the functions Ti, T2, ■■■■Tk are such that 



where 



dQ = ^"MG, 

dQ — col{dTi, dT2, dTx), 
dG — col{dGi, dG2, ...dGx), 



/ 1 

G, 1 

G2 Gi 

G3 G2 






1 







1 



\ Gk Gk-1 Gk-2 
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\ 






Gi 1 J 



(2.7) 



After the integration the system fl2.7p we obtain 

Ti = G'i 



i 2 — Or2 2~ 

^3 = ~ G1G2 + 



G? 



14 — tj4 — OriLT3 + (_r]^(_T2 J 



4 2 



5! ' 3! 



Hence, for the function T,- there are the equivalent representations 



m-j 



j = l,...,K. 



m=l 



where a = (ajm) is some matrix. 

Corollary 2.1 Let us suppose that the functions 



K 



are such that 



K 



^ VjLj =0, Lo = L, 
i=o 

t/ien the constructed system \2.5^) . Ii2.6\} admits the first integral 

K K K 

F{x) = g^^ exp ^ v^T ^ = J] 9?^' exp ^ z/,T„ 
i=i i=i i=i 

where z/q, z/i, z/^ are constants. 

For the particular case when N = 2 and the given curves 

(7,(x) = 0,j = l,2,.,M 
are algebraic curve, then from the Darboux's theory follows that 

^j=Kj{x)gj 



(2.8) 
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thus , the condition on the existence the first integral F takes the form 

M K 

^Q^^jKj + ^UjLj = 0, 
j=i j=i 

For the planar polynomial vector field this condition was deduced in par- 
ticular in [1]. In this paper the following definition is given 
Definition of infinitesimal multiplicity 

Let f = be an invariant algebraic curve of degree n of a polynomial 
vector field X of degree d. We say that 

F = fo + fie+... + /fc_ie^-i G C[x, y, e]/e' 

defines a generalized invariant algebraic curve of order k based on f = if 

fo = f, fi, fk-l 

are polynomials in C [x, y] of degree at most n, and F satisfies the equation 

X{F) = FLp (2.9) 

for some polynomial 

Lp = Lo + Lie + ... + Lk-ie''-^ E C[x,?/,e]/e^ 

which must necessarily be of degree at most d — 1 in x and y. We call Lp the 
cof actor of F. Equations 9]) can be written as 

X(/o) = foLo 
X(/i) = /iLo + /oLi 

< 

^ifk-l) = fk-lLO + fk-2Ll + ... + foLk-1 

The vector field fl2.5p . fl2.6p can be applied to extend the concept of infinites- 
imal multiplicity for the polynomial vector field in M^. 

Proposition 2.1 Let 

g,{x) = 0, x= j = l,2,..,M< N 



8 



are invariant relations of a differential system (S). 
Assume that 

= {91, 92, ...,9m,9m+i---,9n} 7^ 0, 

for arbitrary smooth functions Qm+i, 9m+2, ■■,9n- 
Then the following statement hold: 
System (S) can be written as 112. 5\) . 

Proof. 

Suppose that 

x = X{x) (2.10) 

is a differential system having gi, g2,---,gM as partial integrals. Then 
taking 

<l>j = ^dg,{X), j = l,2,..,M 
Xk = Yd9k{X), k = M + 1,M + 2,..,N 

we get that the system (12. 4p . or, what is the same, (12.51) becomes sys- 
tem (12.101) . Note that in the definition of $j and Xj we have used that 

{91, 92, ■■■,9n} 7^ 0. 

Now we shall study the case when the given number of partial integrals 
isS> N. 

U S = N then the differential system (12.51) takes the form 

= $i{a;^'...,^A/,...,^7v} + .-- + $Jv{^i,--,^A/,..,a;^'} j = l--,N. (2.11) 

Proposition 2.2 The differential system ^2.11\) admits the complementary 
invariant relation 

^^(x) = 0, v = N + l,...,S 

if and only if 



9 



dgi{di) dgi{d2) 
dg2{di) dg2{d2) 



dgiid^) 
dg2{dN) 



$1 

$2 



dgM{di) dgM{d2) 
dgM+i{di) dgM+i{d2) 



dgM+i{dN) $M+i 



dgN{di) dgN{d2) 
dg{di) dg{d2) 

or, what is the same, 



dgnidN) 
dg{dN) 



"^lig--, gM--,gN} + ■■ + ^N{gi--,gM--,g} + ^u{gi--, gM--, gN} = o. (2.13) 

We obtain the proof from the equahty 

dgu{v) = ^u, 

which in view of fl2.ip coincides with (12.131) . 
Below we shall use the following identity 

{/i, /2, /tv-i, gi}{g2, gs, gN, G}+ 
+{fi, f2, In-i, g2}{gi, ga, ■■■,gN, G} + ... ,^ 

/2, In-i, gN}{gu g2, ■■■,gN-i, G}+ 
{/i, /2, /iv-i, G}{gi, g2, ..■,gN-i, gN} = 0. 

The proof follow by considering that f l2.14p is equivalent to the relation 



dgi{di) dgi{d2) ... dgi{dN) 

dg2{di) dgiidi) ... dg2{dN) 

dgN{di) dgN{d2) ■■■ dgN{dN) 

dG{di) dG{d2) ... dG{dN) 



{/i, /2, /s, fN-ii gi} 

{/l; /2, /S; In-I-, ^'2} 

{/l; /2, /s; fN-l, gN} 

{fl, f2, /a? In-i, G} 



(2.15) 



It is easy to show, in view of identity (I2.14p that the Eruguin functions 
determined by the formula 



E 



S+N 

ai,a2, ..,«iv- 



gk}X 



Oi,Q2---ajV-l 



1,2,. .,5 



(2.16) 
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are the solutions of fl2.1ip . where Xai,a2...aN-iy arbitrary continuous func- 
tions: 

$fcL,=o = 0, fc = l,2,....,5. (2.17) 

and 

gs+j = Xj, j = l,2, ..,iV 
The differential system (12. lip in this case takes the form 

S+N 

2'"' ^ ^ {S'ai ) • • •) fl'ajv-i ; -^"^ }-^cn ,a2---ojv-l ('^) (2.18) 

oi,..,ajV-l=l 

where 

-^01,02. ..aiv-i(^) ~ {S'l) 5'2; fi'S) fi'Af}'^al,a2...«JV-l(^)• 
PrOpOSition 2.3 Let gi{x), g2{x), gs{x) S > N are partial integrals of 
a differential system (S). 
Assume that 

= {^1, 92, •••,^7v} 7^ 0, 
then the following statement hold: 
System (S) can be written as 112. 18\) . 

Proof. In fact if we insert (12.160 into (12. lip and considering the identity 
(I2.14P we obtain the require. 

In particular for = 2 we deduce the differential system 

= E^=i ~^j{9j, x} + XsMy, /2 19) 

y = Ef=i ^j{9j, y} + >^s+2{x, y}. 

As usual we denote by M.[x] the ring of all real polynomials in the variables 
X = (xi, X2, , , , , xn)- We consider the polynomial vector field in M^, with 
degree n, i.e., 

V = {v\x),...,v^{x)), v^{x)eR[x], j = l,2,..,N 
n = max{deg{v^{x)), deg{v'^ (x))) 

Definition 2.2 
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We say that {g = 0} C is an invariant algebraic hypersurface of the 
polynomial vector field v of degree n if there exists a polynomial K e M.[x\ 
such that 

dgiy) = K{x)g. 

The polynomial K at the degree at most n — 1 is called the cofactor of 
gix) = 0. 

Definition 2.3 

A nonconstant (multivalued) function is said to be Darboux if it is of the 
form 

f^H\[97{x)l , 

where o"j G C, j = 1, 2, ... 5 arc certain constants. 

Definition 2.4 We shall say that the vector field x — v{x) with invariant 
relations 

gjix) = 0, j = l,...,S>N, (2.20) 

is integrable if it admits N — 1 independent first integrals fi, /2, ..,/jv-i, and 
integrable in the Darboux sense if 

9i, 92, ,9s 

are polynomial functions and fi, /2, ...,fN-i are Darboux functions. 
Darboux proved the following theorem. 
Darboux's theorem 

Si Von connait "^(™'+^)("'+^)----("'+"-~^) — ]\/[^ integreles particulieres algebriques 
de systeme 



dXi dX2 dXn J. ^ r ^ TD r 1 

Li, L2, ...,Ln e R[a;J 



on pourra trouver lemultiplicateur du systeme. 

Si Von connait Mn+r integrales particulieres algebriques du meme systeme, 
on pourra en determiner le multiplicateur et r integrales generales. 

Si Von connait Mn+n—1 — q integrales particulieres algebriques Ui, U2, ...,Uq 
on pourra effectuer Vintegration complete. Les integrales se presenteront sous 
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la forma suivante: 



Ui U2 ■■■Uq 



C2 



Ui U2 ■■■Uq 



a 



n-l- 



In [T^ the following result is proved. 
Jounolous Theorem 

Let V he a polynomial vector field defined in of degree > 0. Then v 
admits ^"'^^^-j]^' + n irreducible invariant algebraic hipersurface if and only if 
V has a rational first integral. 



Proposition 2.4 The vector field \2.11\) with invariant relations l{2.20\) is 
integrable if and only if the vector field the Eruguin functions are such that 



= Hfu f2, ■■jN-i,gj}, k 



1,2,..,^ 



(2.21) 



where A is an arbitrary function. 

Proof. Let us suppose that the vector field v is integrable, then it admits 
the representation [25], 



dfiid,) 
df2{di) 



dfi{d2) 
df2{d2) 



df N-i{di) df N-i{d2] 
di 82 



dfi{dN) 
df2{dN) 

dfN~iidN) 
On 



A {/i, /at-i, *}, 



where A is an arbitrary function. Hence we obtain that 

dgjM = A{/i,...,/Ar-i,5'i} 
on the other hand from (2.2) we obtain that 



By compare both we deduce (12.211) . 

We obtain the reciprocity result as follows. 
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Let us suppose that fl2.2ip holds. Clearly that the condition fl2.14l) holds 
identically in this case. 

By inserting (12.211) in (12.111) we obtain 

v(*) = A( {/i, /a, /s, /tv-i, 9i}{*, 92, 93, 9n}+ 
■■■+ {/i, /2, /a, In-i, 9n}{9i, 92, 93, *})■ 

In view of the identity (12.141) we deduce 

v(*) = {91, 92, 93, ...,9N}Hfi, f2, /s, /jv-i, *} = A {/i, fN-l, *}. 
as a consequence the vector field is integrable. 

Corollary 2.2 Let us suppose that the system C^.ljS is polynomial of degree 
n. 

Then it is Darboux integrable if and only if 



g, \ (2.22) 



X; = k{x)^ J = 1,2,..,^, 



where k, v are arbitrary rational functions and f is a Darboux 's function. 
Proof. 

The proof follows from the fact that in view of (I2.14p . (12.221) we obtain 

dgkM = {k + u) {/, gk}, 
{k + z/){/, gk}\g^^^)=o = 0, A; = 1, 2, .., S. 
Example 2.2. 

We shall suppose that the given invariant relations of the differential 
system (12. lip are the hyperplane 

x^=0, J = 1,2,.., AT. 

We choose the Eruguin functions as follows 

{(^1, ...,(PN-1, X^} 



{(fl,(f2---,<^N} 
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where ^kji^^), ^,(0) = 0, k.j = 1,2,., AT and ifj, j = 1,2,., AT 
arbitrary functions. Hence we obtain that this system takes the form 



are 



We shall study the case when 



j = l,2,..,iV 



^V'i = Ef=i Vkj{x^)dx'', j = 1, 2, .., iV 



{^i,^2---,'^n} 7^ 



(2.23) 



(2.24) 



The differential system fl2:23|) . (12:21) is integrable. 
In fact, by considering that 



k=l 



'k[X'' 



dt if j = N, 
if j ^ ^• 



we deduce the existence of — 1 independents first integrals 

{f^i^)^ElJ-^dx' = c„ j = l,2,..,Ar-l. (2.25) 
It is easy to show that the vector field v(*) in this case is such that 



v(*) = g 



{fi, In-i, *} 



N 



f, 9 = ll'^k{x'). 

J k=i 



For the subcase when the invariant hyperplane are such that 
x'' - am+k = 0, k = l,2,..,N,m = l, 2, M 



and 



M 



iPkjix^) = {xy-\ k,j = l,2,..,N^I/k{x^) = JJ(x*^-a„+fe), 

m=l 

then the first integrals (12.251) in this case take the form 

' fj=^^UtlU'Lli^'-<lm+kf^''-, J = l,..,N 

(Ofc+m)"' 



m+k rT-M" / 

iil=l,l^my"'k+l - ClkA 



(2.26) 
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as a consequence the system fl2.21l) . fl2.22p . 02.241) is Darboux integrable. 
An interesting particular case is the following 



/jq3 r^J 

Hence, by making M — >• +00 we deduce 



= X-' 



n (4 - l)^^'--""^^ ^ sinvra:^!^^!:::!^^. (2.28) 



m=l 



We observe that the differential systems of the type (12.231) appear in 
the theory of the Stackel mechanical system [25], [23]. With respect to this 
system we state the following problem: 

Problem 

Determine the real constants Ki, K2, Kf^_i, L in such a way that the 
hyperplane 

N-l 

x^ = J2 ^j^^ + L (2-29) 
i=i 

is invariant of the system ([MID, ([OlD- 
We solve this problem for the case when 

N = 2, if2 = ^ix^ + y^) + xy, 

2 M 2 

= -A n;ii(^ - 1), ^2{y) = A n(7i - 1) ^'-'^^ 

where a, Oi, om are real constants and A is an arbitrary function. 
The system (12.211) in this case takes the form 

x = X{x + ay) Ylf=ii^ - 1) 

(2.31) 

y = \{ax + y) njli(— -1) 

We require to determine the real values of the constants K and L in such 

a way 

y = Kx + L (2.32) 
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is an invariant straight line of fl2.31l) . 

Clearly that the parameter K must be satisfies the relation 



Hence we obtain that = 1, K2 = — 1 satisfies this relation. For M > 2 
there exist at most four real values of K which satisfy this equation. 

By using the algebraic computer packages we can solve the stated prob- 
lem. In particular for the cubic and quintic system in which 

M 

and 

a = ^5, ai = 1, M = 1 
a = 75, ai = 1, a2 = v^-2, M = 2 

It is easy to show that in this case we obtain respectively 
i^i = 1, Li = K2 = -1, L2 = 

and 

iri = l, Li = 0, K2 = -l, ^2 = 0, 

A3 - -2 - — , -t^3 - -2 + — 

J<". — ^ T . — 3_i_v^ 

-'^4—2 2~' -^4 — -2 "T — ; 

7^3 = -1/2-^, L, = l/2-^ 



2 ) -^0 " 2 ' 
The quintic polynomial system in this case was constructed in [1] 



x = {x + V5y){x^ - l){x^ -iVb- 2)2) 
y = i^x + y){y' - l){y' - {V5 - 2)') 

and admits 14 straight lines. 

The Eruguin functions in this case are: 

^n. = {x + V5y){x^ - l)(x2 - (v^ - 2)2), 
$4+„^ = (TSx + y){y^ - l)(y2 - (v/5 - 2)^), m = 1, 2, 3, 4 

hence the constructed system is not Darboux integrable. 
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For the case when in fl2.3ip 
we obtain the system 



X2 
.... , 



12 (2.33) 

nM fV -I N 

In this case the system admits the following invariant straight line 

x = j, x = -j, j = l,2, ..,M 

y = j, y = -j 

y = x, y = -X 

Clearly, the system (12.331) is Darboux integrable 

If M = +00 then the system fl2.33p takes the form 

i = X n'ili(f^ - 1) = sin TTx 

2 

y = y]17=A^ - 1) = sinvry 

for which the infinity numbers of the straight lines 

= J, a; = -j, j = l,2,..,+oo 
= J, y = -3 j = l,2, ..,+oo 
X + 2m, y = —X + 2m, m G Z 

are its invariant. 

The problem of the determination of the upper bound for the maximum 
number of the invariant straight lines {L{n)) for the polynomial system is an 
open problem. 

It is easy to show that |22] 

2?2 + 1 if n is even, 
2n + 2 if n is odd. 




Lin) > 



There exist the following conjecture 
Conjecture [1] 

L{n) < 3r2 - 1 

This upper bound is reached in particular for n = 2,3, 4, 5. 
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3 Inverse approach for the planar vector fields 

In this section we analyze the Eruguin-Galhuhn theory developed in the 
above section for the case when N = 2. 

The differential system f l2.1ip in this case take the form 

fx = ^i{x,g2} + ^2{9i,x} = P{x,y) ,^ 
\ y = ^i{y,g2} + ^2{guy} = Q{x,y), ^ ' ' 

we set V = (P, Q\ 

For the case when this equations admit the subsidiary invariant curves 
^j(x,y) = 0, j = 3,4, 
The Eruguin functions must be satisfy the relations 

$i{(7,-,(?2} + $2{^7i,(7,} + $,te,^7i} = J = 1,2,. ..5. (3.2) 
Hence the Eruguin functions $m can be determine as follows 

*m = Ef=l \{^){93-- 9m]Vik=\,ki^j9k + (Ss^A^-, gm} + >^S+2{ 9m, y})9, 

9 = nf=i 9j 

(3-3) 

where Ai, A2, As'+2, are arbitrary functions. 

The veracity of this representation we obtain by inserting (13. Sp into (13.20 
and by considering the identity (I2.14p which in this case takes the form 

{9k, 9j}{9i, 9m] + {9k, 9m}{9j, 9i} + {9i, 9k}{9j, 9m} = 0, (3.4) 
Corollary 3.1 

The differential equations can be rewritten as follows 

X = g{x,y)(J2j=i~^ji^^y)^^ + ~^s+2) = P{x,y) 
y = 9{x,y){E-=M^^y)^-~^s-,i) = Q{x,y) ^ ■ ^ 

We deduce the proof by inserting (13.30 into (13. ip by using the identity 

(Ea). 

From these relations we obtain the following proposition 
Proposition 3.1 
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Let gj{x,y) = 0, j = 1,2,. ..,S are the irreducible algebraic curves, then 
the polynomial differential system 6]) admits the Darboux first integral 

s 

/(x,y)=ln(n^7?) 

if and only if in h'j. 3\) 

\ - = uoCj = constants, \s+i = y{f, y}, As+2 = i^{x, f, ], 

where z/q, v are arbitrary rational functions and aj = constants, j = 1, 2, S. 
The proof is easy to obtain from corollary 2.2. 
The system ( 13. 6p in this case takes the form 

X = V{x, F} , . 

_ y = V{y, F} ^'^■^> 

where 

i=i 

We illustrate the above results in the following concrete cases. 
Example 3.1 

In this section we give the results exposed in [25], related with the con- 
struction the planar polynomial vector field with invariant circumferences: 

9i{x, y) = {x- aj)^ + {y - hjf - r| = 0, j = 1, 2, S 

The system f l3.6p under the restrictions 

= = 

takes the form 



i^=- Ej=l - ]lra=l,m^j 9m = P{x, y) 

S S ^ I •> 

V = Ej=i - aj) Um=i,mj^j 9m = Qix, y). 



(3.7) 



or, what is the same, 



y = x{x^ + yy~^ Y.U + + y")''" + - = ^(^' 
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Now we determine the arbitrary functions Ai, A2, •••,As in such a way 
that the above vector field is polynomial of the fixed degree n. 
Corollary 3.2 

Let us suppose that 

A,-GR[a;,y], j = l,2,..,5 
Si = degJ2j=i -^i, 
S2 = max{deg{Yfj=i >^jbj), deg(J2j=i ^j^-j) 

then 

n = max{deg{P), deg{Q)) < 2S — 1 + Si, 

and 

n = max{deg{P), deg{Q)) <2S — 2 + S2, 

if 

From this results we obtain the proof of the following result 
Proposition 3.2 

Every configuration of the circumferences in the plane is realizable by a 
polynomial of the degree at most2S + Si — l or 2S + S2 — 2 where Sj, j = 1,2 
are the degree of the polynomials introduced above. 

In a paper [19] the authors proved that every configuration of cycles on 
the plane is realizable (up to homeomorphism) by a polynomial of the degree 
at most 2(m + r) — 1, where m is the number of cycles and r is the number 
of primary cycles (a cycle C is primary if there are no other cycles contained 
in the bounded region limited by C). 

It is interesting to observe that the upper bound for the degree of the 
constructed vector field is independent from whether its cycles are primary 
or not. 

Now we shall study the case when the circumferences form two nests with 
the centers at the points (0, 0) and (a, 0) respectively, hence 
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f 9j{x,y) =x^ + y'^ -r^ = 0, j = l,2,...,h, 
\ 9j{x,y) = {x-af + -r^ = 0, j = h,...,S 

Clearly, for this case the Eruguin functions are such that 

"^1 = "^-s = = ^i,=^ay Ef=,,+i A, Ut=i,m^2 9m 

$2 = = $ii+2 = = $5 = -4:ay \92 IlLu+i 9m. 

The differential system (13. 7p takes the form 

^ = -2(e5=i \ YiLh+i 9k + eUi+i nL=i 9k)y, ^ 
y = 2(e5=i \ nLh+1 9k + Ef=zi+i h nL 9k)x - 2a T.j=h+i nL 

In particular for the case when 

n = 2/ + 1 = 5 + 1, /i=/, 

-2 Ef=z+i Aj- = a; + 

we obtain the vector field constructed in |25j . 

By designating by Fa{x,y), Fo{x,y) the following polynomials 

( Faix,y) =(x + y-a)U%\iix-a)' + y'-r]) I > 1, 
I Fo{x,y) = Fa{x,y)\a=o. 

we can deduce that the above vector field takes the form: 

X = (^Fo{x, y) - Faix, y)^ y = P(x, y) 

y = - (^Fo{x, y) - y)^ x + aFo(x, y) = Q{x, y). 

This system has the following properties: 

1) has only 3 critical points in the finite plane 

(0,0), (^,0), (a, 0). 

2) the Liapunov quantities a and A for the system are : 
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i) 

r a{0,0) = a{a,0) 
\ A(0,0) = A(a,0). 

ii)) 

A(0, 0) = a' UU ^(«' - ( n;=i(«' - r]) - (-1)' UU 
a\0, 0) - 4A(0, 0) = a' ((2 n;=i {a' - r|) - (-1)' UU ^if " n;=i 



'-(f,o) = anU((i)'--l) 



a2 



,4 



I I I 



^(i.«')=-Yn((f)'-'-j)E n (^-4)- 

j=i 1=1 j=i,j^i 

The circumferences do not intersect if 

so 

A(|,0)<0, 

and, as a consequence the critical point (-, 0) is a saddle. 

It is evident that the other critical points are the stability or non stability 
foci depending on whether k is odd or even. 

Hence we obtain that the constructed polynomial vector field of degree 
n — S + 1 admits S — 21 invariant circumferences. 

The proposition 3.1 we illustrate in the next two examples. 
Example 3.2 

The particular case of the Lienard equation 

X — —h(x)x — ah(x)—h(x) = 
ax ax 



or, what is the same. 



X = y + h{x) 
y = ah{x)—h{x) 

LLJb 
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is Darboux integrable. 

In fact, the first integral F in this case is the following 

F{x,y)=g'[^g^' 
where gi, g2, (Ti, ui are such that 



y)=y+ ^ ^ ^ ^ -hix) 

, , 1 - V4q; + 1 , , , 

giyx, y)^y^ h{x) 



1 - V4q; + 1 



1 + V4a + 1 



2V4q; + 1 ' 2V4q; + 1 

It is easy to show that in this case the Eruguin functions are 

1 + V4« + 1 ^ 1 - V4a + 1 
2v4q; + 1 2v4q; + 1 

Clearly, if 4q; + 1 < then the function F takes the form 



F{x^ y) = {y^ + h{x)y — ah^{x)) exp (^y —a — -arctan 



y + 2^(^) 



It is interesting to observe that if the function h admits the following 
development 

h{x) — X + a2X^ + ... 

then the origin of the given system is a focus. 
Example 3.3 
The differential equation 



z = i {aiQZ + aoiz+ ^ ajkZ^ z'^) 

j+k=3 



(3.8) 



is Darboux integrable, where a — (ajk), j,k — 0,2,3 are real constants 
matrix and 

z = X + iy, z = X — iy 
are the complex coordinate in the plane M^. 
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In fact, the equations fl3.8p are equivalent to the cubic planar system 

X = y{aoi - aio + (ai2 - 021 + 3(ao3 - 030) + («i2 - 021 + 030 - 003) y"^) 
y = x{aoi + flio + (ai2 + ^21 + ao3 + 030) x^ + (012 + 021 - 8(030 + ao3)?/^)- 

Hence, by introducing the correspondent notations we obtain the system 

(3.9) 



X = y{a + bx"^ + cy 



y = x{a + P x^ + •y y'^) . 

We shall analyze the case when c 7^ 0. 
Let Qi, g2 are the functions: 

9jix, y) = i^jix^ - Ao) -y^ + Ai, j = 1, 2 
where Aq, A2, z^i, 1^2 are constants: 

^ 7a — ac ^ ab — P a 

° &7 — c/3 ' ^ 67 — c/3 ' 

Z/^ - 1/2 7^ 0, 



then the following relations hold 

d9jiv) = '2xy{-f-iyjc)gj, j = l,2 
{91, 92} = 4:xy{iyi - V2). 

The proof is easy to obtain after some calculations. 
The given vector field is Darboux integrable with F : 

(z/i(x^-Ao)-7/^ + Ai)^+-^- 

here we use the relation 

c(z/i + z/2) = 7 - &. 

Now we shall study the case when z/i, z/2 are complex numbers. 
By introducing the notations 

7 — 6 = 2cg, 7 + 6 = 2cr 

we obtain that 



(3.10) 
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The system (3.9) takes then the form (we put c — 1) 



( x = y{a+{r-q)x^ + y^) ,^ 
\ y — x{a — {p"^ + q'^) x'^ + {r + q)y'^). ^ ' ' 



By considering that in this case 

gi{x, y) = q{x^ - Aq) - + Ai + ip{x^ - Ao), 
we obtain that the first integral F takes the form 

F{x,y) = (^{y^—Xi—q{x^—Xo)y+p'^{x^—Xoy^exp(^2rarctg 



p{x^ - Ao) 



y2 - Ai - q(x'^ - Ao)' 

4 Planar differential system with one invari- 
ant algebraic curve 

In this section, by applying the results of the section 3, we construct the 
analytic planar vector field 

x = P{x,y) 
y = Q{x,y), 

where P and Q are analytic functions on the region G C M^, from a given 
set of trajectories: 

9jix, y)^y- yj(x) = 0, j ^l,2,...,S >2 

{gi, 92} ^ y'2{x) - y[{x) 0, (4.2) 
nf=i%(^)7^0 

where ^ = y'j and yi, y2, . . . ,ys e C'iG C R), r > 1. 

We shall study the particular case when yi, 1/2, ...yg are solutions of the 
equation 

g(x,yj)^0, j^l,2,..,S 
where g{x, y) = is an algebraic irreducible curve. 
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By considering (13.11) we obtain that the require system fl4.ip . fl4.2p : 

y = y)y2ix) - $2(x, y)y[{x) 



(4.3) 



The condition (13.21) . on the existence of the complementary partial inte- 
grals, in this case take the form 

<^i{y2{x)-yj{x))y+<!>2{yj{x)-yi{x)y+<l>j{yi{x)-y2{x)y = 0, j = 3,4 . . . , ^. 

(4.4) 

The given set of differential equations (13. 6p in this case can be rewritten as 
follows 



a; = Ef= A n^^jd/ - ym) + ~^s+2g = P{x, y) 

y = Ef=i hy'j Ylm^j^y - ym) - ~^s+ig = Q{x, y) 

or, what is the same. 



X = y'Xs+2 + y'-\EU \ - E;=i 

y'-\Y,U - ^1 vi ^lA) + ••) + ••• 
y = -y''^Xs+i+ z/^"HEf=i ^jy'j ^^s+i Ef=i yj)+ 

, y^~\J2j=i \y'j - J2j=i yj Ef=i Xkyky'k-) + •••• 



(4.5) 



(4.6) 



we set V = (P, Q). 

Proposition 4.1 Let us suppose that the arbitrary functions \i, ...A5, A5+1, A5+2 
are such that 

= -go (a;), A5+2 = Po{x) 
s 



= ^C£,Pk{x)y'f\x)) 



Ao 
A, 



fc=0 

s 



(4.7) 



A;=0 



where pk, Qk o,re continuous functions on D G 

1 1 ... 1 

Ao 



yi y2 



yj 



yr 



yr 



yr 



1 

ys 



yf 
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A, 



1 1 

yi y2 



yr yf 



1 1 ; 1 

yj-i yj+i '■■ ys 



5-2 S-2 ■ 5-2 



Then the differential system takes the form 

X = Poix)y^ + ....+ psix) = P{x, y) 
y= qo{x)y^ + .... + qs{x) = Q{x,y) 



(4.8) 



and we set v = (P, Q). 

We shall study the case when po. Pi, ■■■,Ps, Qo, Qi, ■■■Qs are polynomials 
on the variable x and such that v represented a polynomial vector field of 
degree 

n = max{degP, degQ) 

Corollary 4.1 Let g be a irreducible polynomial on the variables x and y: 

s s 
9 = aoix)Yl{y-yj{x)) = ^a,(a;)/-^ (4.9) 



j=0 



where 



ai = -ao(x) X;f=i% (a;) 
a2 = ao{x) Uj<kyji^)ykix) 



(4.10) 



^ as = i-l)^ao{x)llLiyji 



[X 



then 



dg{v) = {J2^)9 = K{x)g 
j=i 9j 

where K is the cofactor oi g = 0. 

Proposition 4.2 Let the curve Iji4-^li4- ^ ) ^■5 invariant of the non zero 
polynomial system of degree n : 



X = ri(x)y" ^ + .... + r„(x) = P{x, y) 
y = go(a;)2/" + .... + Qnix) = Q{x, y) 



(4.11) 
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where rj{x), qj{x), j = 0, 1, .., n are polynomials of degree j in the variable 
X, then 

S < 2n. 

Proof, al absurd, let us suppose that S = 2n + l, then from (4.9), (4.7), (4.11) 
we obtain that 

Pn+l{x) = 0, 

Pj{x) = 0, 
?i(a;)=0, j = 0,l,..n, 

on the other hands from (14.71) . (14. 8 p we obtain that 
~Xj{x) = 0, J = 1,2, ...,2n + 2 
hence the vector field is a zero vector field. Contradiction. 

5 Quadratic system with one invariant alge- 
braic curve 

In this section we shall study the case in which the vector field (14.51) is 
quadratic i.e., 

X = Ps-2y^ + Ps-i{x)y + Ps{x) = P{x, y) 
y = qs-2y'^ + qs-iix)y + gsix) = Q{x,y) 

where max {degP, degQ) = 2 and qs-j, Ps-j, j = 0, 1, 2 are polynomials 
in the variable x. Below, for simplicity we shall denote this system as follows 



i= Poy'^ +Piix)y + p2ix) = P{x,y) ,^ 

y= qoy'^ + qi{x)y + q2{x) = Q{x,y) 

First we prove the following general results related with the system (15. ip . 

Proposition 5.1 Let us suppose that / I5. 1\) is such that 

dg{v) = {aoy + ax + (3)g, (5.2) 

where g is given in the formula (I4.9l) . (l4.10p and ao, P are real constants 
and Pj, qj are polynomials of degree j in the variable x : 

Pj ~ Ylk=oPjkX^ ^ 3^ 
Qj=Ei=o1JkX^ j = 0,1,2. ^ • ^ 

29 



Then, 



If Po 7^ hence 



max{deg aj{x)) <j, 
ma.x{deg g) < S. 



If 



hence 



Po = 0, pu 7^ 0, 
ao = {Sk + m)pn, go = kpn 
Sk + meN 



max(deg aj(x)) < kj + m, j = 1, 2, .., 
m8ix{deg g) < Sk + m.) 

In fact, from (15.21) we deduced the differential system 

A - — = B ■ Si 

dx 
dao 
Po-r- = 

d 

P2—^ + q2as-i = {ax + l3)as 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



dx 



where 



a = col{ao, ai, as) 
is a vector and A, B are matrix which we determine respectively as follow 







f Pi 


Po 













\ 










P2 


Pi 


Po 
























P2 


Pi 


Po 






























P2 


Pi 


Po 










\^ 











0P2 


Pi, 


J 
































ax + P — Sqi 




-is 


-1) 


qo 

















-Sq2 


ax 


+p- 


-{s 


- 1 


)qi 


tto - 


(s- 


2)go 



















-2g2 ax + l3 - qi ao 
-q2 ax + (3 j 
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From fl5.7p we easily deduce that if po 7^ then the coefficients aj, j 
0, 1, 2, .., S' are polynomials of degree at most j. 

For the second case, after integration we easily deduce that 



Rp 



(5.8) 



where 



m+fc 
1 Pi ; 



P = coI{pT,pT^\ 
is a vector and R is the following constant matrix 



m+kS\ 
Pi ) 










































Rl 


R\ 


































Rl 


Rl 













Rl 


d2 

-"-fc+l 


















Rl 


Rl 


i?2 









; 


R3 d3 
^k-'^k+l 


^k+2 


00 


^2k 













JXq 


Rf 








Rl 


R-k 


T^S 

^k+1 


^k+2 











R 



Hence we easily obtain the veracity of our assertion. 

Corollary 5.1 // the invariant algebraic curve is non reducible and k = 
> then m = 0. 

Pn 

In fact, if m 7^ and under the indicated condition we have that the given 
invariant curve is reducible. 

Proposition 5.2 The maximum degree of the irreducible invariant algebraic 
curve of the non Darboux integrable quadratic system ( 15. ij) is 12. if pu ^ 



In fact from (15. Sp . (15. 5p by considering that the given curve is irreducible, 
then we obtain that m = 0. 

On the other hand from the last of equation of (15.71) in particular we 
deduce that if /c > 3 then 



— ^lfe(PioP22 -P10P21P11 +P20Pn) = 0. 
Pu 



(5.9) 



Hence, if = then the degree of the algebraic curve is not maximal. 
On the other hand if 

pIqP22 - P10P21P1I + P20PU = 

then the quadratic system is Darboux integrable. 



■Sk 
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Corollary 5.2 Let us suppose that the algebraic curve 

S kS 

g{x,y) ^J2J2^'M^''+P^^y^"y' = (^-10) 

1=0 j=0 

is invariant curve of the maximum degree of the quadratic vector field. 
Then: 

v= {{piiX+pio)y+p22x'^+P2iX+p2o)d^+{3puy'^+{quX+qio)y+q22x'^+q2iX+q2o)dy 
and 

12 

dg{v) = — {pliy + P22X + P21P11 - 2p22Pio)5'- 

This results can be extended analogously for the polynomial system of degree 
n. 

Proposition 5.3 Let us give the invariant curve with S branches of non- 
Darboux integrable polynomial system of degree n : 

with cof actor K = ^^=0 where aj{x) = ^i=oTjkx'^, and 
T — (Tjk) is a constant real matrix. 
Then 

maxde^(^) < { I) , ^ [J' , ^ (5.12) 

yyy^ ~ ^ S{n + 1), tfpo = 0,pn^O ^ ' 

For the case when pn = it is easy to show that 

"0 = Sqo, 
max deg Uj < j , 
max deg g < S. 

Prom this result and proposition 4.2 we deduce the proof of the following 
proposition. 

Proposition 5.4 The maximum degree of the invariant curve of the non 
Darboux integrable polynomial planar vector field of degree n (^{n) ) is 2n{n-\- 
l), i.e., 

N(n) < 2n(n+ 1). (5.13) 
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6 Quadratic system with one invariant alge- 
braic curve. Examples 

First we study the case when S = 2, i.e., we analyze the quadratic system 
with invariant algebraic curve of the type 

g{x, y) = ao{x)y'^ + ai{x)y + a2{x) = (6.1) 

where Oq, ai, 02 are polynomial on the variable x. 

We shall study the cases when po = 0, pn 7^ and po ^ 0. 
For the first case the he equation (15.101) takes the form 

2 2k 

g{x,y) = 5^^R;(piix + pio)V = 0. (6.2) 
1=0 j=0 

Clearly, the maximum degree of this curve is six. 

We shall illustrate this particular case in concrete examples. 

Example 6.1. The quadratic vector field 

x = ax'^ + {-2ac + y + b)x + 2hCia^ - lOaCirf + C^d^ 
y = 3y2 + {-Viae + 66 + dx)y + (-ISOCia^Cg - Sa^ + 60CiadC3+ 
da - 6CiC3rf2)x2 + [dh - 2adc)x - TSOa^C^CsS + 225Cia^ + AbOa'^C^dC^i 
-2d^Ci - imCiO^d + 12a2c2 + Sfe^ - 90^2^20^3 - 12ca6 + 29Cid'^a + Qd^C^Cl 

where Ci,C3 are arbitrary nonzero constants, admits the invariant curve of 
degree six 

1 3 

g = Ciy"^ + {x^ + {Cid - 3Cia)x - ACiac + 2Cih)y + x^C^ + {-d - -a)x^+ 

(6 - 2ac)x3 + (-aCsC^d^ _ 9^^^^ ^ 21Cia2 - TSCgC^a^ + Cid'^+ 

3QC^Clad)x^ + (-3Cia6 - 2Cidac + Cidh + Qa'^Cic)x+ 

7^ 1 c; 

l50Cs3C^a'^d a^C^d + —d^C^a+ 

i 1 2^2^ 

Cilb^ - 250C3Cfa^ - -d^Cl + 2C3Cld^+ 

12s ^ 

a^Cf - ACiahc - ZQC^Cld^a + Aa^c^Ci = 

2 

with cof actor 6(?/ + ax — 2ac + 6). 
Example 6.2 
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Now we study the case in which S = 2, po 0. Clearly, in this case the 
invariant algebraic curve of the quadratic system is the family of the conies 

g = a^y^ + (anx + aio)y + 022^^ + 0210; + 020 = (6.3) 

where ao, an, aio, 022, 021, O20 are real constants. 
In particular, for the quadratic system 

i; = _ ^ 2)a;t/ + (/? - 4)x2 - 2/?x 
ij = {a- A)y'^ - 2{a + 2)xy + (a - 4)x2 - 2a x 



(6.4) 



we have that po = (3 ^ 

After some calculations we can prove that the invariant curve is the 
parabola 

{y -xf -2x = Q 

with the cofactor 

K = 2{(3 + a)y + 2{(3 + a)x - 2(3. 

The quadratic system with invariant parabola was constructed in [Sadl] 
and admits the equivalent representation : 



i = P {{y ~ ^)'^ — 2x) — Ax {x + y) 
y = a {{y — xY — 2 x) — 4 (x + y) 



2 Af^< ..\2 (6-5) 



The points 

0(0,0), M,|,^) 

are its critical points, where = 2[[a — f3y — 2a). The bifurcation analy- 
sis show that this system is generic. The bifurcation curves divide the plane 
{a, (3) in 17 region in which we observe a qualitative change in the behavior 
of the trajectories of the constructed quadratic system. We determine 38 dif- 
ferent quadratic systems, among these there is one with one limit cycles [21] . 
The quadratic system 

X = + By^ + 2C)Ao + j{-2Byqu - 2xBq22 + xA^A2Bq2i)y 

y = ^(^1^^ + ^^11^ + '?22a;^ + ^2ia; + AiC), 
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admits as invariant the curve 

Ax^ + By^ + 2C = 

with cofactor K = Aqx + Aiy, where A, B,C, Aq, Ai,qu,q2i,q22, are real 
constants such that A ^ 0, B 0. 

Now we study the case when S = 3, i.e., we analyze the quadratic system 
with invariant algebraic curve of the type 

2 3A; 

g{x,y) = XlE^^K^^ii^ +^^io)V = 0. (6.6) 

1=0 j=0 

Example 6.3 (The Filipstov system ). 
For the quadratic system 

r X = 16(1 + a)x - 6(2 + a)x2 + (2 + 12x)?/ 

\ y = 3a(l + a)x^ + (15(1 + a) - 2(9 + 5a)x)y + 16?/^ 

we have that 3k = 4 =^ m = 0. 

As we can observe in this case the quadratic system possesses the irre- 
ducible invariant algebraic curve 

1 3 3 

g{x, y)=y^ + -(3(1 + a) - 6(1 + a)x)y'^ + + + + a)a'^x^ = 0. 

The cofactor of this curve is 

K = A8y- 4(1 + a)x + 5(1 + a). 

Example 6.4 For the quadratic system 

f x = {2 + 3a)x^ + (2 + 4y)x + y 
j y = (5 + 4(1 + a))y + 6?/2 + ax'^ 

we obtain that k = 2 =^ m = 0. 

The quadratic system possesses the irreducible invariant algebraic curve 
of degree four [6] 

a^x"^ + +2ax'^{x + l)y + (1 + x)y'^ + y^ = 
The cofactor K in this case is 

K = 18y + {5 + 6a)x + 5. 
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For the case when S = A we obtain the invariant algebraic curve of the 
type 



2 4A: 

g{x, y) = J2Yl ^^5(^11^ + p^<^yy' = 0- (6.7) 

1=0 j=0 

Clearly, the maximum degree of this curve is 12. The upper bound is reached 
in particular in the following example [6j 

Example 6.5 For the non Darboux integrable quadratic system 



X = xy + X + 1 
. ^2 81 2 57 

?/ = 3?/ - yX + — 

has we have that /c = 3^m = 0asa consequence the above vector field 
has the invariant irreducible algebraic curve of the maximum degree 12. In 
the indicated paper was showed that the the curve 

-442368 - 7246584a;2 + 71546517x^ - 97906500a;'' + 41343750x« - 234375003;^°+ 
48828125x^2 ^ (322272x - 12126312x3 + 23463000x5 + 1125000x^ + 15625000x9)?/- 
(98784 - 711288x2 + 5058000x^ - 375000x^)?/2 + (32928x - 1124000x3)?/3 - 5488?/^ = 

is its invariant. 



7 Construction the polynomial planar sys- 
tem with invariant algebraic curves with 
variables separable 

In this section we deal with the polynomial system with invariant algebraic 
curve with variables separable 

(7(x,y) = Fi(x) + F2(y) = 0, 
where Fi, F2 are arbitrary polynomials : 

deg{g{x,y)) = max(^deg{Fi{x)), max{deg{F2{y))J . 
We state and study the following problem. 
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Problem 7.1 Let g he a function: 



Y[{x - aj)dx + B Yliy - ^j)dy 

where A, B, ai, 02, ...,am,i, &2, ai'e real parameters such that 

oi < 02 < .... < ami, ^1 < ^2 < bm2, AB ^ 0. 

We require to determine the non-Darboux integrable polynomial vector v of 
degree n for which the given curve is its invariant. 

We propose the solution of the state problem for the following particular 

cases 

772,2 = 71 — 1, rrti = n — 1, 
m2 = 2m — 1, mi = 2m + 1, n = 2m + 1, 
m2 = 2m — 2, mi = 2m, n = 2m, 
m2 = mi = m, n = 2m + 1, 
m2 = mi = m, n = 2m + 2. 



m,2 



Proposition 7.1 T/ie polynomial system of degree n fEj 

( x = {Ax + By + C)dyg{x,y) = P{x,y) 

1 y = -{Ax + By + C)d^g{x, y) + \g{x, y) = Q{x, y) 

admits as invariant curve 

/n—l „n— 1 

- hj)dy + K2 - aj)dx = 0, 

i=i i=i 

where ai, 02, ...,a„_i, 61, b2,---,bn-i, Ki, K2, go, A, B, C, A are arbitrary real 
parameters. 

By choosing the arbitrary parameters properly we can construct the non- 
singular algebraic curve of degree n, hence the genus (G) of this curve is: 

G = i(n-l)(n-2) 
Example 7. 1 Let g is a. nonsingular curve of degree 2m + 2 such that 

/m „ m 

Y[x{i-f - l)dx + l[y((y-f-l)dy = 
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It is easy to show that the polynomial system of degree n — 2m + 1 : 

X = {A^x + B^y + Cm) n7=i y{{—f - 1) = y) 

y = -{AmX + Bmy + Cm) n^i x{{—)'^ - 1) + >^9m{x, y) = Q{x, y) 

admits as invariant the given curve. 
By considering that 

oo 

TT((— )^ - 1) = sini/ ^ lim 51^(0;, y) = 5(0 + cos a; + cost/ 

and choose the arbitrary parameters Ami ^m, Cm properly, we obtain the 
analytic planar vector field 

X = siny = P{x, y) 
y = — sin X + X{gQ + cos x + cos y) = Q{x, y) 

for which the curve 

go + cos X + cos y = 

is its invariant. 

Clearly, the constructed analytic system admits infinity many number of 
limit cycles. 

Proposition 7.2 The polynomial vector field of degree n 

r x = {a + byx)dyH{x, y) ,^ 

\ y = -{a + byx)d^H{x,y) + {n + l)byH{x,y) 

admits as invariant the algebraic curve of degree n + 1 

H{x, y) = x^+' + G'„_i(x, y) = 0, (7.2) 
where is an arbitrary polynomial of degree n — 1. 

Clearly, this system in general has no Darboux integrating factors or first 
integrals. 

The following particular case is an interesting one: 
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Corollary 7.1 Let 

px m+1 „y m—1 

9m{x, y)^go+ Yl ^(^^ ~ + / n y^y^ ~ ^f^'^y ^ ^ 

Jxo Jyo 

is a curve of degree 2m+2 with the maximum genus G = 2(to + l)(m — 1) + 1 
is invariant of the vector field of degree n — 2m + 2 : 

x = {a + brnyx) njL~/ yiy"^ - b]) 
y^-{a + b^yx) HJlV ^(^^ ~ «i + (2^ + '^)bmygm{x, y) 

Example 7.2 By making m — > +oo and choose the arbitrary parameters 
properly we deduce from the above system as a particular case the analytic 
system 

r x = aJo{y) 

\ y = -aJo{x) + Xy{Ji{x) + Ji{y) + go) 

where Jq, J\ are the Bessel functions. This analytic system admits an 
infinity many number of limit cycles. 

Analogously we construct the polynomial system of degree n — 2m : 

x= {a + b^yx)Y[JZ,\y^-b]) 

y^ - (a + bmyx) YY^^i {x^ - a^) + 2mbmygm{x, y) 

with invariant curve 

„^ m+1 „y m,— 1 

gm{x,y)=go+ JJ (a;^ - a|)rfa; + / Yl{y^ - b'^j)dy ^ 

Jxo j^i Jyo j=i 

is a curve of degree 2m + 2 with maximum genus G = 2m{m — 1) 

By using the algebraic packages it is possible to show the following propo- 
sition. 

Proposition 7.3 There exist polynomials p{x, y) , q{x, y) of degree n for which 
the non-Darhoux integrable differential system 

r i = adyg{x, y) + p{x, y), a = const. ,^ 
1 y ^ -adxg{x,y) + q{x,y) 
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has the invariant curve of degree n + 1 



g{x, y)^g, + A ([[{x - aj))dx + B ( J](y - bj))dy = 0, (7.4) 




for certain values of the real parameters go, A, B, ai, a^, bi, 62, b^, 0,0, bo. 

Clearly if this curve is non singular then the genus is G = ■::n{n — 1). 
Example 7.3 

Let us suppose that the given algebraic curve (7.4) is such that 



It is possible to construct the non-Darboux integrable polynomial vector 
fields of degree n = 2m + 1. In particular for n = 3, 5, 7 we construct the 
following non-Darboux integrable polynomial systems. 

For the polynomial system of degree seven 



my\^ - 32p2g6 + 32y2g6 ^ iQOpYq^ + 6Ay^q^)iyo + \y{-12p^x^ - Ap^y^ - Ay^q'^+ 
Sp^y'^x'^ + Ap^y"^ + — p^q^ + I2x^p^q^ + Sy'^x'^q^ — Sp^y'^q^ + Sx^ + 2x'^p'^+ 
2a;2g4 - 87/^^2 - 12x\^ - pY + + ^^yV) 
y = -l/6Ax{-6ApY - 32p2g6 _ qq^Y + 32p^x^ + 32x'^q^ - 32p^q^- 
192xVg^ + 160j9^a;2g2 + leo/^^g^ + 54^2^6 _ gg^4^4 ^ 64a;6g2)/(p2 ^ ^2^^^ 

— l/64a;A(8p^|/^a;^ — 12/;^?/'^ — Sy'^x'^ — p^q"^ — %x^p^q^ + Sy'^x'^q^ + 12p^y'^q^ + 2p*y'^— 
4a;V - '^x^ff + +P^ + '2y^q'^ - p'^q^ + '^P^x'^ + Ax'^q^ - I2y^q^ + Sy^) /(p^ + q^) 




the invariant curve is 



g{x,y) = 1/8X^ + (-l/4p2 _ l/4g2)x6 + (l/Sp^ + l/2p2g2 ^ 1/8^4)3,4 _ ^(^4^2 _^ ^2^4^)^2^ 
l/Sy* + (-l/4p2 _ l/4g2)y6 + (l/8p4 ^ l/2p2^2 ^ l/8g4)^4 ^ (_l/4p4^2_ 

l/4p2g4)^2 _ i/i28g8 + l/32gV + 13/64gV + l/32gV _ i/i28p8 = 



with the cofactor 



K — \{—y + x){x + y){—p^ — + y"^ + x^)yx. 
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In this example we have that 

For the quintic vector field 

\ X — — i/o(6r^ — 6y^)(r — 3y^)y — A(r^ — x^)x(r^ — 3y^) 
The curve 

111111 1 

/ \ 9z] 4.9 fi 4.9 4.9 

q(x,y) = -v r -y H — r u H — x x r H — r x r = 

^^'^^ 6 3 6 6 3 6 39 

is its invariant with the cofactor 

X = 6y2A(3x2-r2). 
The cubic polynomial system 

r x^ uo{y^ -y) + \x{2q^x'^ - 2p2y2 _ 3^4 ^ ^4-) 
\ y = -i/o(a:^^ - a:) + Xy{2q'^x^ - 2p^y^ + 3g^ - p^), 
has as invariant curve of degree four 

g{x,y) = l/Sq* + 1/8/ + l/4a;^ - l/2gV + 1/4|/^ - l/2pV = 
with cofactor 

K = 4A(2g2a;2 - 2p2^/2 _ ^4 ^ ^4^,^ 

It is interesting to observe that the above constructed polynomial system 
of degree seven under the change 

[ X = VX, y = ^/Y, 

[ ^ ~ I y ' X )\-=^^y=^^ 
can be transformed to the cubic system 

Y = (-64p2 - Uq^)X^ + (96g^ + 192^2^2 ^ 96^4)^2+ 
(-32p6 - 32g6 - 160gV - 160gV)X + 64^^?^ + 32p2g6 + 32^6^2)^3 ^ (.gys^ 

(I2p2 + 12g2)y2 ^ (_2p4 _ 2g4 _ gp2^ _ ^2j92^2 _ 8^^2 ^ 8X2)F + gV _ ^6 _ ^6^ 

X = -((-64p2 - 64g2)y3 + ^9g^4 ^ i92p2^2 ^ 95^4)^2^ 
(-32p6 - 32g6 _ l60gV - 160gV)>' + 64^4?^ + 32p2g6 + 32^6^2)^3 + 8X^+ 
(12p2 + 12g2)x2 + (-12p2g2 _ 8p2y _ §^^2 ^ §^2 _ 2^^ _ 2p4)X- 
p^-q^ + gV + - 4y2g2 + 4yg4 + 8p2rg2 ^ 4^4^ _ 4p2y2-)^) 
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which admits the invariant curve of degree four 

g{X, Y) = 1/8X^ + (- 1/V - + (l/2p2g2 + 1/3^4 ^ i/Sq^)X'^+ 

(-l/4gV _ i/4gy)x + 1/8Y^ + (-l/4p2 _ l/Aq^)Y^ + {l/2pY+ 
l/8p^ + l/8g^)y2 + (-l/4gV - l/4gV)F - l/128p^ + l/32p^q^- 
l/128gS + 13/64pV + l/32p2g6 

with cofactor 

K(X, Y)^(X-Y)(X + Y- (p' + q'))X. 

Example 7.4 

Finally we analyze the case when 

9{x, y) = ^(2^0, yo)+ n ^(^^ ~ ^ / n ^(^^ ~ 



is a curve of degree 2m + 4. 

Analogously to the above case we can construct a non-Darboux integrable 
polynomial vector fields of degree n = 2m + 2 

In particular for n = 4, 6 we construct the following polynomial systems. 

For n — 4: 

X = (1/6|/V - 5/108|/2 - l/18a;^ + 7/324x2 + l/729)i/o + (1/6^ - l/54:xy)X 
y = (-l/54xy + l/4:xy^ - l/12yx^)uo + (l/4y^ - 5/108x^ - 1/18^^ + 4/2187 + 4/243x2) A 

the invariant curve is 

g{x, y) = l/Qx^ - 5/36x^ + 2/81a;2 + l/4y^ - IjlSy"^ + 4/2187 = 0. 

The cofactor in this case is 

K = i — (1458x*^ - 1215a;^ + 216a;2 + 2187|/^ - 486^^ + 16) A. 

236196^ ^ ^ 

For n = 6 we construct the following vector field 

' X = (-l/128yV - 4/3y2 + 1/64x6 - l/2x^ + l/12y^ + 40/9^2 + l/8yV- 

256/27)z/o+ (l/8xy3 - 3/512x 

• / 1 5 5 3 1 5 

y — I xy H xy + —yx — 

^ / 96 18 ^ 48^ 

Q?/a; + — x2/)z/o + (-— 1/ 

11 4 1 4 ro24 2 2 1 

+ 4^ ^ ^ " 2^ + + I^^ ^ 
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The invariant curve and its cofactor are respectively 

g{x,y) = l/67/'^-16/37/^+128/32/2+l/8a;^-16/3a;6+704/9a;^-4096/9a;2+65536/81 = 0, 



K = {-Ay + y^ - 3/6Ay^)X + {1/8 - 8/3 x^ - l/16y*x + y^x + 32/3a;)z/o. 

We observe that after the change 

x = X^ y = Y\ 
vix) wix) 



V 

' x y 
the above system takes the form 



X = Y{{-1/2X^Y^ + l/SX^ye _ 3/512X^Y^^)iyo+ 
[-1/128^^X^ + 1/64X12 - 256/27-4/3^4 + 40/9X4 + 1/8^4X4- 1/2X«+ l/12r«)A) 
Y = X((-l/128Fi2 - 2Y'^ + 1/4^^ + 1/16X^2 - 11/6X8 + 16X4 - 1024/27)z/o+ 
(-l/96X2yio - 4/9^^X6 + 5/18X2y6 + i/48y2xio + i6/27X2y2)A) 

This polynomial system of degree 13 admit as invariant curve of degree 

16 



g{X,Y) = l/6ri2-16/3r8+128/3r4+i/8Xi6-16/3Xi2^704/9XS-4096/9X4+65536/81 = 
with cofactor 



K = — (i/o(24Xi°-512X6+192r4x2-12rSx2+2048X2)+A(192F6-9ri°-768F2)rX) 
96 

We observe that the constructed above curve are singular curves. 
We state the following problem 

Problcma To construct the differential system (7.3) for which the in- 
variant curve (7.4) is non singular. 

Clearly, if the curve (7.4) is non singular then the genus is 

G = -n(n-l). 
2 ^ ^ 
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as a consequence the maximum number of algebraic limit cycles are -n{n — 
1) + 1. 

To construct the require vector field should be satisfied the relation dg(y) = 
K[x)g under the condition that the curve g{x,y) = is non singular. To 
obtain the explicit expression for the vector field in general it is necessary 
to solve a lot of technical and theoretical problems. In particular for the 
cubic system we have 31 parameters which must be satisfies 27 equations. 
By solving these equations we obtain that the nonsingular curve of genus 3: 

' g{x, y) = l/Ax^ - l/3(a2 + ai)x^ + l/2a2aix^ + 1/%^ - 1/3(62 + bi)y^ + l/262&iZ/'+ 

l/12/i(^ - hhlhlal - hhlhlal + 26^62^2 + 26264^2 _ 2hia\a% + 2hihlal+ 
^ 2blbfal — 2hyji\aia2 + 2afalbl + 5blblaia2 — 2b2b\aia2+ 

hbia\alb2 + 2blala\ - Safa^fo^ - 2bia\alb2 + 2ala\bl - hbla\al + 2bla\ai^ = 0, 

= (^{aj - aia2 + ai)(-&2&i + b^ + bl)j 

which is invariant of the cubic system if the parameters ai, 02, &i, 62 are roots 
of the certain homogenous polynomials P, (ai, 02, &i, 62), j = 1, 2, 3 such that 

f Pj{\ai, \a2, Xbi, A62) = X^Pjidi, 02, bi, 62), j = 1, 2 
\ -PslAai, Aa2, A61, A62) = X^^Psia^ 02, &i, ^2) 

8 The Poincare problem and 16th Hilbert prob- 
lem for algebraic limit cycles 

8.1 The Poincare problem. 

The question on the existence an effective procedure to find a natural 
number N(n) which bounds the degree of all irreducible invariant curve of 
a non-Darboux integrable polynomial system of a degree n is well known as 
Poincare problem [20] . [3] . [25] . 

The problem on the existence of the polynomial planar systems with an 
invariant algebraic curve of the maximum degree was studied in particular 
in [5], 

citeSad2. In [5] the author gave an explicit polynomial system of degree n 
for each non-singular real algebraic curve (? = of degree n what is that 
system's invariant. The author states also that n is optimal for a generic 
class of algebraic curve. 
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Proposition 8.1 Let N(n) is the maximum degree of the irreducible alge- 
braic curve (4-9, (4-10) that is invariant curve of the polynomial system (3.5) 
of degree n. 
Then 

{n + l)< N(n) < 2n{n+l) 

The proof of the lower bound follows from the proposition 7.2 and 7.3 and 
the upper bound from the proposition 5.4. The upper bound is reached in 
particular for n = 2. 

9 The 16th Hilbert Problem for Algebraic 
Limit Cycles 

In 1900 Hilbert [13] proposed in the second part of his 16th problem to esti- 
mate a uniform upper bound for the number of limit cycles of all polynomial 
vector fields of a given degree. This question has been studied in particular 
in [2Bj,[26] for algebraic limit cycles. 

By considering that the ovals of the invariant algebraic curve are isolated 
periodical solutions of the vector field for which is it invariant we deduce that 
the maximum number of algebraic limit cycles of the polynomial system of 
degree n with one invariant curve (we denote this number by 1) ) is at 
most the genus of the curve +1, i.e., 

A{n,l) < G + 1 

hence we observe that to solve the 16th Hilbert's Problem for Algebraic Limit 
Cycles it is necessary firstly to solve the Poincare problem for this case, i.e., it 
is necessary to find the maximum degree of the invariant curve [25] . 

Proposition 9.1 

A{n, 1) < (2n^ + 2n - l){n'^ + - 1) + 1. 

The proof follows from the proposition (5.4) and from the Harnak theorem 
on the maximum numbers of the ovals of the algebraic curve. 
From the above results we prove the following 
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Proposition 9.2 The maximum number of the algebraic limit cycles for the 
polynomial planar vector field of degree n with one invariant algebraic curve 
A{n, 1) is such that 

A{n, 1) > max{hQ{n), hi{n), h2{n), h^^n), h^i^n), ) 

where hj{n) are the maximum numbers of ovals of the algebraic curves 
H,{x,y) = 0: 



^m„,n—m 

1 



Hi{x,y) = aa;"+i + E^io "n-™™^"?/""". 
Hkix,y) = JiUT'^ix - aj)dx + Uiiv - bj)dy), 
ai < a2.... < flmfe, hi <h2.... <hi^, A; = 2, 3, 4 
m2 = I2 = n — 1, = n + 1, I3 = n — 1, 7714 = n, 1^ = n. 

From this inequality we deduce the following result |26j : 

Corollary 9.1 

A{n,l) > -(n-l)(n-2) + l. 

Conjecture 

. / -, X nin — 1) 
A{n, 1) > ^ ^ ' 

This conjecture can be solve if we construct the polynomial system of degree 

n with one invariant nonsingular irreducible invariant curve of degree n + 1 

(see proposition 8.1). 

It is interesting to observe that in [16] stated the following problems: 
Is 1 the maximum number of algebraic limit cycles that a quadratic system 

can have?. 

Is 2 the maximum number of algebraic limit cycles that a cubic system 
can have?. 

The answer to this questions for the system (3.5) with S > 1 was given 
in [26]. 

Is there a uniform bound for the number of algebraic limit cycles that a 
polynomial vector field of degree n could have?. 

What is the maximum degree of an algebraic limit cycle of a quadratic 
polynomial vector field?. 
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What is the maximum degree of an algebraic hmit cycle of a cubic poly- 
nomial vector field? 

The partial solutions of these problems we obtain from the above results 
and results proposed in [25], [26] . 
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